REPORT 956 


LINEARIZED COMPRESSIBLE-FLOW THEORY FOR SONIC FLIGHT SPEEDS 
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SUMMARY 

The partial differential equation for the perturbation velocity 
potential is examined for free-stream Mach numbers dose to 
and equal to one. It is found that, under the assumptions of 
linearized theory, solutions can be found consistent with the 
theory for lifting-surface problems both in stationary three- 
dimensional flow and in unsteady two-dimensional flow. 
Several examples are solved including a three-dimensional 
swept-back wing and a two-dimensional harmonically oscil- 
lating wing, both for a free-stream Mach number equal to one. 
Momentum relations for the evaluation of wave and vortex 
drag are also discussed. 

INTRODUCTION 

Much of the recent progress in the theoretical analysis of 
compressible-flow fields is attributable to the successful 
application of linearization methods. Although the basic 
assumptions used in conventional linearized theory appear 
at first glance to be highly restrictive, it has been found that, 
just as in the analogous case of thin-airfoil theory for incom- 
pressible flow, the methods have many fields of utilization 
adequate for most engineering purposes. Since the basic 
methods are so well known and depend on such relatively 
simple mathematical tools, it appears obvious that the 
range of applicability of the theory should be explored. 
Such is the purpose of the present report. 'It has been more 
or less tacitly presumed in the past that such applications 
cannot treat cases for which the flight velocity is near the 
speed of sound. In the study of two-dimensional steady- 
state problems in airfoil theory, this presumption is certainly 
true. The Prandtl-GIauert and Ackeret rules for variation 
of pressure coefficient with free-stream Mach number in the 
subsonic and supersonic regimes, respectively, are clearly 
invalid for Mach numbers near one, since perturbation 
velocities become arbitrarily large. In this case, linearized 
theory therefore predicts its inability to treat such problems. 
On the other hand, if linear methods are applied to non- 
stationary two-dimensional airfoil and particular steady- 
state, three-dimensional, lifting-surface problems at sonic 
speeds, a consistent theory results since solutions are found 
which yield perturbation velocities of the same order of 
magnitude as those calculated for free-stream Mach numbers 
of, say, 0.6 or 1.5. 

Unfortunately, arbitrary thickness distributions at sonic 
speeds cannot be studied by linear theory in the steady state 


since, in general, the theory predicts infinite pressure differ- 
ences between the wing surface and infinity. In the particu- 
lar case of a yawed, symmetrical wing of infinite aspect 
ratio, the results are, however, again consistent with the 
theory and yield pressure distributions which are the same 
as those determined by using only the component of free- 
stream velocity normal to the leading edge. The derivation 
of this latter result for a free-stream Mach number of one 
will be given. 

The difficulty of not being able to include thickness 
effects in general, together with the uncertainty of the mag- 
nitude of the viscous effects, leaves the question as to the 
limitations of such a linear theory in application to practical 
wing shapes. Such a question can certainly not be resolved 
by mathematical reasoning alone. The extent to which the 
fluid medium can be idealized at these speeds is left, for 
the time being, unsettled and it remains for experiment to 
determine whether the consistent mathematical results ob- 
tained from the linearized equations provide reasonably 
exact predictions. In this connection, it should be mentioned 
that the few experimental results available for the total lift 
on thin triangular wings at Mach numbers near one tend to 
confirm the theory. But even if more detailed experimental 
results indicate that further refinements are necessary, there 
is still little doubt but that the linear potential solutions will 
provide a valuable basis for more exact extensions of theory. 

The present report is divided into three parts. In the 
first part, the linearization of the partial differential equation 
for the velocity potential is carried out in some detail for 
steady-state conditions. A byproduct of this derivation is 
the nonlinear form of the equation for two-dimensional flow 
used by von Karm&n (reference 1) to determine his similarity 
rules for transonic flow. The equation for unsteady two- 
dimensional flow based on the same assumptions is also given. 
The second part of the report is restricted to two-dimensional 
unsteady problems for values of Mach number near one. 
The principal contribution of this section is the evaluation 
of the change with time of the pressure distribution over an 
airfoil starting suddenly from rest at a speed close to that 
of sound. Such an idealized problem involves a step func- 
tion in velocity in which the airfoil has zero velocity for all 
negative values and near sonic velocity for all positive values 
of time. From these results the initial build-up of lift can be 
calculated from Mach numbers near one, although the 
eventual value of the lift cannot be found by linear methods. 
Further application can also be made to problems in flutter 
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and gust loads. The third part of the report treats the 
steady-state three-dimensional problem. Both lifting sur- 
faces and symmetrical nonlifting wings are considered and 
it is seen that in the former case consistent solutions are 
obtained by particularly simple means. . ' These solutions rep- 
resent the limiting case of both subsonic and supersonic lift- 
ing-surface theory and give, for example, the same value of 
lift-curve slope at the speed of sound . that was obtained for 
the supersonic triangular wing by Stewart (reference 2). 
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LIST OF IMPORTANT SYMBOLS 

y coordinate of trailing edge, y=a 1 (x) 
free-stream speed of sound 
local speed of sound 
critical speed of sound 

aspect ratio [~ . ( s P an ) — "j 
r Lwmg areaj 

Fresnel's cosine integral cos ~x 3 dxj 

wing root chord 

lift coefficient f , . ^ r-1 

Lowing area)J 

do, ,, 

da 



indicial lift coefficient 

Cb. 

vorte* drag coaBceat 

D 

drag 

D . 

vortex drag 

erf(x) 

eiTor function of x e d\j 

E(fn,kn) 

elliptic integral of second kind 


(J 0 Vl— £n 2 sinV<i>) 

E n 

*(%*■) 

F( in,K) 

elliptic integral of firet kind 


/ CK d<(> V 


vjo Vl~^tt 2 sinV/ 

K n 

'(y*.) 

K«',E n ' 

elliptic integrals with moduli k' n 

K 

modulus of elliptic functions 

k'n 

"" ill 

l 

length of body 

Mo 

free-stream Mach number (— ^ 

Ni (at, v) 

„ ■ . ■ - - - • • 

I [cos x s C(ut—z 2 ) — sin x 2 S (cut — a: 2 ) ] dx 

J V' 

IV* (cut, v) 

f [cos x 2 S(ait — x 2 ) +sin x 2 C(ut—x 2 )] dx 

*/ yy - .. 

Ap 

loading coefficient (pressure on lower sur- 


face minus pressure on upper surface 
divided by^fre e-stream dynamic 
pressure) . 

free-stream dynamic pressure ^ p o y o s ^ 


S 

s 

So 

V 
t 

to 

u,v,w 

V 
F 0 

x,y,z — - 

a 

p 

y 

A<f>o,Au 0 ,Aw e 


Pa 


4 

cu' 


u 


polar coordinates in yz plane (y~r cos 9, 
g=r sin 6) 

Fresnel’s sine integral (X sto f *’*) 

operational equivalent of t 

wing semispan 

time 

a 0 t' 

maximum distance measured parallel to y 
axis from x axis to trailing edge (fig. 15) 
perturbation velocity components in x,y,z 
directions, respectively 
local velocity 
free-stream velocity 

radial component of perturbation velocity 
Cartesian coordinates 
angle of attack in radians 

Vji=W[ 

ratio of specific heats, for air y=\A 
discontinuity in component in 2=0 plane 
semivertex angle of swept-back wing 

CxJCq 

2 

free-stream density 
total velocity potential 
perturbation velocity potential 
impressed frequency (reference to truo 
time) 


co_ 

a a 


I— THE LINEARIZED EQUATIONS OF MOTION 

STEADY STATE 

The nonlinear partial differential equation satisfied by the 
velocity potential <£ of an isentropic flow field can be expressed 
in the form 


2 2 2 


( 1 ) 


where the subscript notation is used to indicate differentiation 
"and a is the local speed of sound given by the relation 


(C) = Mo ! { 




m-'j 


( 2 ) 


In this latter equation V 0 and M 0 are, respectively, velocity 
and Mach number of the free stream, y is the ratio of specific 
heats (for air, 7= 1.4), and V is local velocity. 

Introducing the perturbation velocity potential <f>, where 


0= — V<jX+$ 


(3) 


it is possible to express equation (1) in terms of the deriva- 
tives of 4> and the parameters M 0 and V 0 - To begin the 
linearization of the resulting equation, the coefficients of the 
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second ordered derivatives of $ are expanded in Maelaurin 
series with ascending powers of 

■^r, —■ The convergence is assured provided 

fO n h “ 

<1 (4) 

or, in a slightly modified form, provided 

[ V s — F 0 2 | <~?j == 5a 0 I (5) 

U V / U) 

If the assumption is now made that T^Cl so that 

Vo Vo Vo 

second and higher powers in the perturbation velocities can 
be neglected in comparison with one, the partial differential 
equation can be simplified to the form 

0xx j 1 — M 0 2 [~L (y — 1)M 0 2 | + <£,,*+ 4>t*— 

24 > xt £ M*- 2 4 >*y A/.*- o (6) 

Vo v Q 

From this equation all the succeeding expressions will be 
derived. 

Two- and three-dimensional linear equations, M 0 ?*l. — 
Since equation (6) is obviously nonlinear, additional assump- 
tions must be made to reduce it to a linear form. Clearly, 
these assumptions must involve the relative magnitudes 
of all the terms in order to determine which ones may be 
neglected. Perhaps one of the least restrictive set of con- 
ditions is that: 

(a) The ratios of the perturbation velocities to the free- 
strearn velocity are small enough to be neglected when 
compared to one. 

(b) The nondimensional velocity gradients are also small 
in comparison with one where, in the nondimensionality, 
perturbation velocities are measured in terms of free-stream 
velocity and the length is measured in terms of a char- 
acteristic length of the body. 

With the aid of these assumptions, it follows that, to the 
order of the approximations made, the perturbation velocity 
potential <j> satisfies the well-known linear equation 

(1 — Mo 2 ) 'frzx J r'f>vv J r4>u=Q (7) 

In the case of two-dimensional flow, the equation is inde- 
pendent of y and thus may be written in the form 


[Or— 1) , ii 2 +r 2 +w 2 \ 

| 2 VFo + Fo 2 ) 


( 1 — Mo 2 ) 4>xx-{-<i>zi=0 


( 8 ) 


Two- and three-dimensional nonlinear equations, M a = 1. — 
The study of equation (8) in both subsonic and supersonic 
flow has shown that for arbitrary lifting surfaces or sym- 
metrical nonlifting airfoils the value of the induced velocity 
u on the surface of a fixed geometric configuration is pro- 
portional to (|l— M 0 2 j) _1/2 . In all airfoil problems, the value 


of u becomes infini tely large as M 0 approaches one, either 
from above or below, and the basic assumptions are thus 
violated. Such a difficulty led Oswafcitsch and Wieghardt 
(reference 3) and Sauer (reference 4) to abandon the restric- 
tion of linearity and to seek a more exact equation at M 0 =L 
Retaining the assumptions underlying equation (6) and 
setting F 0 =a* where a* is the critical speed of sound, it 
follows that at M 0 =l the perturbation velocity potential 
satisfies the equation 

4>x4>xx — <i>zx -r 4>x4>xx— 0 (9) 

Since <f> x is much larger than <f>, as the Mach number ap- 
proaches one, equation (9) may be further simplified to 

<t>x4>xr~4>xx= o (10) 


If, in three dimensions, the perturbation velocities do not 
remain small, equation (6) again supplies the necessary 
form of the differential equation at M 0 =L From the rela- 
tion F 0 =a*, the required expression is 


(r+1) 

a* 


4>X<f>IX 4>vr <£z 2 +„* 


2 

<M>*v=0 


( 11 ) 


Two- and three-dimensional linear equations, M a —1 . — 
Equation (10) has been used by von K&rm&n (reference 1) 
to establish similarity rules for two-dimensional transonic 
flow and is the basis for work continuing at the present time. 
(See also reference 5.) If at M 0 =l the assumptions made 
in the linearization process still hold, it follows from either 
equation (10) or (8) that the differential equation reduces to 
the form 

4>«=0 (12) 

It is possible, however, to predict independently from this 
relation that linearized methods cannot be applied to the 
calculation of arbitrary airfoil pressure distributions. The 
range of applicability of such an equation is thus almost 
nonexistent. On the other hand, the linearized form of 
equation (11) or (7) at M 0 =l is 

= Q (13) 

and from this equation a class of nontrivial solutions can be 
obtained for particular boundary conditions. Both equa- 
tions are of parabolic form in the number of dimensions for 
which they are defined. In the present report, formal solu- 
tions satisfying the imposed conditions will be obtained in 
three dimensions for flat lifting surfaces with swept-back 
leading edges and for an infinitely long, symmetrical, swept- 
back wing. 

UNSTEADY STATE 

The derivation of the steady-state equations for the veloc- 
ity potential was developed in some detail because of the 
various results to be obtained. Similar methods can be used 
when unsteady conditions are to be considered, the diiferen- 
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tial equation for the velocity potential being now in the 
form 


— ($ <' f 4 * 2 $* < ' + 2<b v /' + 2 $„ t >) + 

2 2 2 

^2 ^2 $»*$*$* ^ssy^*^|/ = 0 


(14) 


where t' represents time. The details of the derivation can, 
however, be avoided by referring directly to the equation 
satisfied by the velocity potential for the propagation of 
sound waves of small amplitude. (See reference 6, p. 492.) 
In this form of the equation the Cartesian coordinate sys- 
tem x, y, z is assumed fixed in the medium so that free-stream 
velocity is zero, while the wing, which moves in the direction 
of the negative x axis with velocity Vo, generates small pres- 
sure disturbances. As a consequence, the velocity potential 
of the field satisfies the well-known wave equation in three 
space dimensions: _... - ... . 


q 2 ‘Pt'f 4*xx iii 4>zt — 0 ( 15 ) 


Equation (15) is reducible to canonical form by means of 
the relation 


t=a 0 t' 


and the three-dimensional form of the equation is therefore 

4>tt <i>xx <i> w $n— 0 (16) 

while in the two-dimensional case independence with respect 
to y yields 

4>tt $zx $11— 0 (17) 

II— TWO-DIMENSIONAL LINEAR PROBLEMS FOR 

Mq near one 

UNSTEADY STATE. Aft 5 1 

It was pointed out in the derivation of equation (12) that 
the linear equation for the velocity potential is not applicable 
to airfoil problems in either the subsonic or supersonic 
regimes for M 0 near one. The possibility still remains, 
however, of analyzing unsteady flows during the period in 
which the perturbation velocities remain small. As an ex- 
ample of such a problem, consider the case of a flat lifting 
surface at a small angle of attack a starting from rest at 
a velocity V 0 near the speed of sound. The perturbation 
potential for such a motion is equivalent to the change in 
potential brought about by an abrupt change a in angle of 
attack of an airfoil flying in a steady-state condition at ve- 
locity equal to V 0 . The change in load distribution brought 
about by this maneuver is the so-called indicial load distribu- 
tion and a knowledge of such indicial functions is important 
for applications using operational calculus. 

Figures 1 (a) and 1 (b) furnish an insight into the nature 
of the boundary conditions for airfoils traveling, respectively, 
at supersonic and subsonic speeds. The chord is initially 



Ca) Supersonic wing. 

(b) Subsonic wing. 

Figure X.— Boundary conditions for two-dimensional unstcndy-lift problem. 


on the x axis with leading edge at the origin and trailing 
edge at jr=c B . With increasing time, the wing section 
travels in the negative x direction and sweeps out a portion 
of the xt plane as shown in the figures (indicated by shaded 
areas). Throughout this part of the plane the boundary 
conditions require tjaat the induced vertical velocity is 
— V 0 a, while elsewhere on the plane the induced velocities 
are continuous functions of z. In figure 2, sketches of the 
airfoil in the supersonic and subsonic cases are shown 
together with indications of the manner in which the dis- 
turbance field spreads. These wings are presented in zyz 
space with time as a parameter so that their coordinate 
system is not to be confused with that of figure 1. The 
airfoils are traveling from right to left at Mach numbers of 
0.8 in. the subsonic case and 1.2 in the supersonic case, and 
for a time corresponding to that required for the wing to 
travel a distance of one-third chord length. At t = 0 cylin- 
drical waves are induced at each disturbance point, that 
is, at each point of the chord. These waves expand radially 
at the velocity of sound but, because of the velocity of the 
wing, the center of the expanding waves moves relative 
to the initial disturbance point. At a given instant in time 
the entire disturbance region of the airfoil is contained 
within the closed surfaces shown in the figures, the outer 
surfaces corresponding to the largest values of time. In the 
supersonic, case, the pressure distribution over the wing 
reaches a steady-state value as soon as the wing moves 
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FlGUTti 2. — Sketch showing extent of disturbance fields after travel of one -third chord length. 


Region B (between lines x=—t, z=t, and x=c Q —M a t) 
A P- f - arc cos . V+W~l ( ! ± arc sin HYI 

<Z - VW=I It arC C0S x+AU 1 tM 0 \2 ‘ arc sin t)\ 

(18b) 


Region C (between lines x=t, t— 0, and x=c 0 —M 0 t) 


A p _ Act 


(18c) 


From the pressure distributions it is possible to calculate 
the indicia! lift coefficient C La (f) as a function of Mo and t. 
Since 


C,At) 


_1 ft* 

Co Jo 


A p 

qa 


dx 


ahead of the expanding cylindrical wave produced at t = 0 
by the leading edge. In the subsonic case, the wing never 
leaves the disturbance field of the cylinders and, as will be 
seen later, the steady-state pressure distribution is 
approached asymptotically. 

It is apparent from equation (17) that the characteristic 
cones have semivertex angles equal to 45° and that the 
cones with vertices on the xt plane have traces with slopes 
equal to ±1. These cones determine the upstream bound- 
ary of the field of influence of the vertex point and their 
cross sections in the plane t = constant are the disturbance 
regions of the cylindrical waves arising at the vertex. Thus, 
perturbations in pressure produced initially at the leading 
edge of the wing section are confined at later time, in the 
modified coordinate system, to the cone with vertex at the 
origin and traces x= ±f. 

The solution of equation (17) for boundary values of the 
type under discussion has been indicated in reference 7 
through consideration of an analogue problem in supersonic 
lifting-surface theory. Thus, the shaded areas in figures 1 
(a) and 1 (b) are thought of as swept-forward lifting surfaces 
situated in a stream directed along the positive t axis at a 
Mach number M 0 = y2. The boundary values remain the 
same; that is, <t> I =w= — V B a on the wing and 4>t, 4>v, 4>t are 
continuous functions of z elsewhere in the xt plane. In 
lifting-surface terminology, the unsteady case for supersonic 
speed becomes a wing with supersonic leading edge, while 
the case indicated in figure 1 (b) involves a subsonic leading 
edge. 

The solution for the wing traveling at supersonic speed 
has been given in reference 7 in a form valid for all Mach 
numbers greater than or equal to one. The expressions for 
load coefficient A plq, where 

A P = Pl~ P* 

2 5P0F0 2 


differ analytically in various regions of the xt plane. These 
expressions are: 


Region A (between lines a ;=— M 0 f, x=—t) 


A p 4a 

T - VM 0 2 -i 


(18a) 


the following results are obtained: 


First time interval 0<f< ^ ^°M ~ 




C Cn 

Second time interval W <£< 0 


1-f-Mo x x Mo — 1 
. c a —Af 0 t\ 


(19a) 


^®=7[i(2+ are9 “‘ JL -r E )+ 

1 £+M 0 Co — tAI B , 1 — t — — ryrx*"l 

7M?n arc cos - — * +se 


(19b) 


Third time interval 


Co 

Mo-1 


<f 


C,(f) = 


VMo*— 1 


(19c) 


These results have been discussed in reference 7 for values 
of M 0 greater than one. They still hold, however, for sonic 
flight speeds and, in fact, can be reduced to the expressions: 


First time interval 0<£< ^ 

C L (£)=4 (20a) 

cc 

/» 

Second time interval <£ 

C ‘.®-Kl +“ c ^“ 5! r+ 2 V^) < 20b > 

The indicial lift coefficient is seen to be constant and equal 
to 4a up to the time t' =^~ or up to the time required to 

travel one-half chord. Following this first time interval the 
indicial function rises monotonically, reaching an infini tely 
large value as time increases. The growth of CzJS) is, of 
course, in agreement with the fact that the steady-state load, 
coefficient becomes infinitely large in linear theory for M 0 =l. 


« 
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This means that the theory cannot be used to predict the 
complete extent of the C L Jf) variation with time but that 
during the earlier part of the motion the assumptions remain 
valid. In figure 3, curves of C L Jt) are plotted as functions 

SM 0 t 



Figure 3.— Indiciai-lift-curve slope for Mach numbers between 0 and 1.4 shown to time re- 
quired to travel 12 half-chord lengths. 


of for. values of M 0 — 1 as given by equations (20a), 

(20b), and M 0 =1.2, M 0 =1.4 as given’ by equations (19a), 
(19b), and (19c). Also included in the figure are variations 
of C L Jf) for M 0 =0 as calculated from Wagner’s results 
(reference 8) by R. T. Jones (reference 9) and also for 
M o =0.8. The derivation oj results leading to the Mo=0.8 
curve will be given subsequently in this paper. The value 
of C!t a (f) at M 0 =0.4 for a short interval of time is also 
drawn. The dashed portions of the curves were not calcu- 
lated but were drawn to agree with the known asymptotic 
value of the lift function. . 

Application of the indicial lift function at M 0 =l, — Once 
the indicial lift function is known, it. is possible to determine 
the lift corresponding to a given variable motion. Consider, 
as an idealized example, the case where the airfoil experiences 
an abrupt rising and sinking motion atjregular time intervals. 
Such a motion involves abrupt plus and minus angles of 
attack without rotation or pitching so that a(t') is given by 
the meander or square-wave function shown in figure 4 (a). 
In this example the variation of a is such that the curve for 
C h {t r ) can be calculated easily. In figure 4 (b), C L {t') is 
shown for the case in which the discontinuities occur at 
intervals of time equal to Cq/Vq, that is, after each chord 
length of travel. The principal point of interest in this 
example is the fact that such a motion yields no excessive 
value of lift or perturbation velocities and the entire analysis 
is within the frame-work of linear methods. 

When the variable motion is more complex in character the 
lift coefficient can be expressed by means of Duhamel’s 
integral. Corresponding to the angle-of-attaek variation 




Cb) Resulting variation of lift. 

Figure 4.— Lift resulting from square-wave angla-of-attaefc variation. 

a(i') as a function of time, the lift coefficient C L (t r ) is given 
by the expression 

C L {i/)=^^C La {t'-r') a {r')dr' (21) 


In analysis related to equation (21) it is convenient to employ 
techniques associated with the use of the Laplace transforma- 
tion. (See reference 10.) Thus, if the Laplace transform 
/ (s) of the function fit) is defined by the relation 

/(«)=J^ e- ,l fit)dt 

then equation (21) can be rewritten in the form 

^is)=sC r Js)a{s) ( 22 ) 

Consider now the case of a lifting fiat plate oscillating har- 
monically without pitching at a frequency « and maximum 
angle of attack equal to a ma!C . Setting 

«(t)=a max e 1 “‘ = a max e < “ a »" (23) 


then equation (22) yields 


cm 

®max 





(24) 


By straightforward manipulation, the inverse transformation 
of equation (24) can be shown to give 

+ e -ty [duf-v) -i S(,i-r)l- 

2-y/^ i A 7 ! {cat, ?)— 2-y^ A/j [cot, v) | — — arc tan c °° 

(25) 
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where 


OIC , j ' co'Co 

v ~ 2 2 a 0 

os' is true impressed frequency [a— a maz e ia ' 1 ’) . 

C{ast—v), S(ost—v) are Fresnel’s integrals (reference 11). 

iVi(&>f, v), Ntfat, v) are integrals defined in the list of symbols. 



Figure 5.— Lilt resulting from cosine-wave angle-of-attack variation. 


If the response to a cosine variation of ct is required, only 
the real part of equation (25) is used. Such a response in 
the early stages of the maneuver is shown in figure 5. For 
very large values of time C L as given by equation (25) 
approaches the value 

£rr 4e< "‘ (Vi ( 26 ) 

from which both the amplitude and phase shift of 0 L resulting 
from either an impressed sine or cosine variation of a can be 



Figure 6.— Amplitude o t oscillatory lift resulting from a cosine angle-of-attack oscillation 
(without pitching) at.\fc=l. 

readily determined. Figure 6 shows the amplitude of lift 
oscillation corresponding to continuous angle-of-attaek oscil- 
lation, as determined from equation (26), plotted as a func- 
tion, of v=~^- It is apparent that as os approaches zero, 



C L approaches infinity; that is, as the impressed wave ap- 
proaches the “step” function (fig. 7), the lift coefficient ap- 
proaches infinity. This result is in agreement with equation 
(20b) as t approaches infinity. As the frequency parameter 
v is increased, however, the value of CV la „ is reduced and 
reaches a minimum of about 3.4 for a value of r=0.9. For 
a speed of sound around 1,000 feet per second and a wdng 
chord of 6 feet, this would correspond to a frequency of 47.7 
cycles per second, a value well within the range of practical 
flutter frequencies. It is interesting to note that figure 6 
also shows that as the frequency of oscillation becomes large 
(i. e., r>3) the value of mjir approaches the value 4. 

This is the same as the value for C L Jt) in the early stages 
following a step variation of a. 

UNSTEADY STATE. Aft<.l 

It has been pointed out that, in the determination of the 
indicial lift function for a wing traveling at subsonic speeds, 
the lifting-surface analogue involves the calculation of load 
distribution over a swept-forward wing with subsonic edges. 
This means that recourse cannot be made in the solution to 
the simple source distribution method used in reference 7 
to treat the U/ 0 =£ 1 case. Since, however, a portion of the 
leading edge in the present case is still supersonic, the prob- 
lem is particularly adapted to lifting-surface methods de- 
veloped by Evvard in reference 12. Figure 8 indicates, as 
in figure 1 (b), the geometry associated with the boundary 
conditions. In the Eward analysis the solutions, as in the 
previous case, are calculated for various regions. As an aid 
in identifying the different results the sketch denotes these 
regions by Roman numerals. 

The load distributions and indicial lift functions will be 
given here for values of t up to the time when the character- 
istic from the trailing edge first crosses the leading-edge 
trace, that is, for 


or 


c 0 






Co 


do (1 — A/”o) 
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c 0 

T=W 


Figure 8.— Regions used in the study of subsonic unsteady lift. 

This period covers the time in which the wing travels 
2Mo/(l-M 0 ) half chords and, since the present analysis is 
concerned with values of Mach number near one, will in 
some cases extend beyond the range of the linear theory. 
The following results are obtained for load coefficient: 

Region I (between lines x—t, t— 0, and x—c a — t) 

Ap 4a- . 

qW a (27a > 

Region II (between lines (x— — M 0 t, x=t, and a:=cb— f) 

A p_ 8 a ( M 0 / t—x , _ £ lM 0 t+x\ , 

2 ~^M 0 \1+M 0 V Mot+r +arc tan V ~t=TJ (27b) 


Region III ^between lines x=c 0 —t, x=t, x=-t-{- -^^- r 
and 

1-Mo) 

Ap 8a / t—x i 4a f . 2x—t(l—M 0 ) , 

2 ~V(1+M 0 ) V x+M4 + AMo L arc sin (1+M„)t + 

Region IV (between lines x=t, a?— c 0 — M 0 t, and x=c Q —t) 

. Ix-Cp+AU m 

2 ^Jl 0 V t(M 0 -l) (27d> 

Region V ^between lines x=c 0 —M 0 t, and 

t [f 

2a . x , 4a . 2x — t(l— M 0 ) . 32 aK 

— r -r arc sm T+— it arc sm — H — 
ttM 0 t irM„ (1+M 0 )£ 1 x a (I+Mo) 

I c a —x—Mtit 2 a . 2c 0 — f(l+M 0 ) r Ar 

V(T tii+M-3 +if ‘ 


where 


k'=JT=F 


k ~yl l (t+x) 2 (l+Mo) 

, . IX + M) 

ip=&vc sini - 

M Co 


N = P ^ r 1 - f . _j (tt+cj (1— M 0 )[(t?— iii) (1+M 0 )— c 0 V2]+2g 1 [g(l+M^ —Cq)2—u{1—Mq)] \ , 

3 friMoJo LV(«-yi)(M+«i)J ( [» (1 + M 0 ) — c 0 V2 + Cj (1 —M 0 ) ] [« (1— M 0 ) —Ui ( 1 +M 0 ) ] ) 1 


where . ... 

«-*) 

v= 7f (i+x) ’ ■■■ T' ■ ' 

k) j i QCOm pIeie elliptic integrals 

K, E complete elhptic integrals 

In figure 9 (a) the growth of pressure distribution with 
time is shown at subsonic speed for the period of time 
covered by equations (27). For purposes of comparison, 
pressure changes calculated from equations (18) are shown 
in figure 9 (b) for supersonic flight velocities. 

Equations (27) suffice to determine the initial growth of 
indicial lift coefficient at subsonic speeds. Such results 
were given in figure 3 at M 0 ~0.8 along with the calculated 



(a) Subsonic, 
(b) Supersonic. 


Figure 9. — Pressure distributions on wings receiving sudden anglc-of-attacic change at 1—0. 
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growth for about one chord length of travel at M a =0A. The 
value of C La (t) at f= 0 is, for all flight speeds, equal to 4/M a . 

Expressions for C L J$) are as follows: 


First time interval 


Co 

1 + A/q 




(28a) 


Second time interval i —M a 


C ‘« (4) i m arc tan V 2 f(l+M a )- 2 C( + 

V ( f -lj|ro) + ( 1 l+M^ V[2i(l + M 0 )-2 Co ][2c 0 -^l-A/o 2 )H- 


4c 0 — 2f(14- Mg) 

i+M, 


arc sin 


“ V i(1 <a+1 h ‘ v^(i+j«-v+ 


[2o»-i(l+AA)] arc tan V tH+^-c S+l/V^. (®Dr* 

1+Af, 


(28b) 


where (^^^is given by equation (27e). 

HI — THREE-DIMENSIONAL LINEAR PROBLEMS FOR M t 
NEAR ONE 


STEADY STATE 

General solutions for arbitrary Mach numbers. — Two 
methods of attack are available for the solution of linearized 
problems at sonic speeds. In the first place, solutions to 
equation (13) can be written formally and the extent to which 
these solutions satisfy the original assumptions can then be 
investigated. In the second place, general solutions of equa- 
tion (7) can be studied in the limits as M 0 approaches 1. 
Since this latter method furnishes added information con- 
cerning the variation of the variables with M 0 , it will be used 
first. 

In linearized theory the boundary-value problems of wing 
theory are concerned with two separate properties of the 
wing; the thickness effects and the effects produced by the 
twist, camber, and angle of attack. The first is called the 
nonlifting case and the second is the lifting case. Solutions 
of equation (7) for 1 are given in reference (13) as 
follows: 

In the nonlifting case 





Aw e (x u yQ dxidyi 

V(*-*i y-Pto-yd'-P* 


(29) 


where p—-yJM a 2 —l and Av) 0 =2w 0 where w a is the vertical 
perturbation velocity on the wing and therefore related 
directly to the slope of the wing surface relative to the x axis. 
The integration region t is the area on the wing within the 
Mach forecone from the point x, y, z. 

In the lifting case 


<£(x, 2 /, 2 )=-^ 

where A<f> 0 is the jump in the value of the velocity potential 
in. the plane of the wing. The sign 1 denotes “finite part” 


1 / 


dxdyi 

[{x-xtf-fny-ytf-FzY 1 * 


(30) 


of the integral and introduces special integration techniques. 
(See reference 13.) 

Equation (29) expresses the velocity potential for the sym- 
metrical wing in terms of an integral involving supersonic 
source distributions while equation (30) employs doublet 
distributions. In the two cases the distributions are de- 
termined from the geometry and the load distribution over 
the wing, respectively. 

Source and doublet distribution effectiveness at infinity. — 
It is well known that the lift, drag, and pitching moment of a 
given wring may be calculated either from direct integration 
of the local pressures on the wing or by means of momentum 
considerations where the induced velocities of the wing are 
determined at an infinite distance and the desired forces are 
related to an integration over a control surface enclosing the 
wing. In the three following sections the latter approach 
will be considered and the limiting value of drag at M 0 =l 
computed. The initial portion of this theory requires the 
evaluation of source and doublet effectiveness at infinity 
and the concept of equivalent source position, an idea which 
appears to have been given first by W. D. Hayes in reference 
14. 

Consider, as in figure 10, a point P with coordinates x, y, z 
lying within the induced field of a supersonic wing. The 



Figure io.— C oordinates used in study of supersonic source. 
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Mach forecone from P is given by the relation 
Z— Zi=|8V(2/— y0 2 + ( 2 — Si) 2 


(31) 


where x u y u Z\ are running coordinates of a point on the 
surface of the cone. Introducing polar coordinates 

y=r cos 6, z=r sin d 

and rewriting the abscissa of P in the form 

x=x 0 +fir 

it follows that the trace of the forecone in the zi=0 plane is, 
in the limit as r approaches infinity, 


%\=&yi cos 6+x o 


(32) 


It is, moreover, possible to show that the effect on the veloc- 
ity potential at the point P as r approaches infinity is the 
same for all points (x u y u 0) for which x y —^y x cos 0= constant. 
The value of this effect is 


<j>= 


2x^2pr(x 0 —x 1 +py i cos 9) 


(33) 


and follows from the asymptotic evaluation of the supersonic 
source potential 


<t> z 


1 


2iryj(x—x l ) i —^(y—y 1 ) i —^z i 


for large values of r. The potential at P for the source at 
(a?i, y u 0) is thus the same as for the source shifted along the 
trace to —Pyi cos 6, 0, 0), the intercept of the trace on the 
x axis. For Mach numbers near one, equation (33) can be 
rewritten 


4 >= 


2ir^2^r(Xo—Xi) 


(34) 


and is equivalent to the potential at P for the source at 
(®i, 0, 0). The induced velocities at P due to a source at 
(a?i, y\, 0) follow immediately, for arbitrary M a and for M 0 
near one, from the gradients of <t> in equations (33) and (34). 
It is important to note that equation (33) is a function of 
the azimuthal angle of P so that, in general, a source does 
not have a fixed equivalent position with respect to its po- 
tential at infinity; equation (34), however, is independent of 
the azimuth 6. 

The source-sink potential is applicable to the study of 
symmetrical nonlifting wings.. "When lifting surfaces are to 
be analyzed, the doublet potential 


0= 


_ _ fi 2 g - - 

~2t[(x— Xi) s — 0 1 {y— 3 /i) a — £ 2 s 2 l s/i 


must be considered and the question of equivalent doublet 
position with respect to the potential at infinity arises. In 
this case the doublet position can again be shifted parallel 
to the trace of the Mach cone from P at infini ty and the 
potential at P is given by the expression 


j8 2 s 


~2ir[2l3r(xo—Zi+f}yi cos 0)] 3/2 


(35) 


and, for Mach numbers near one, 

£ 2 2 


4>= 


27r[2j3r(x 0 — 2 i)] 3/j 


(30) 


Momentum relations. — The vectorial force F on an aero- 
dynamic body inside a control surface S is given by the 
surface integral 

” J* J ^ ^ ~ IJ f>V [( f)'ds] (37) 

where vector notation is used and 

0 subscript indicating free-stream condition 

p,p local static pressure and density 

V local perturbation velocity vector 

For the purposes of the present report, equation (37) will 
be modified according to the assumptions of linearized theory 
and the surface S restricted to a semi-infinite circular cylinder 
of radius r, its axis of symmetry lying along the x axis, and 
with one face in. the x=0 plane while the other face is at 
x — constant . (See fig. 11.) 

From the linearized theory, 


u 


and 


Po 1 0 


p-p 0 =-pu |Vo«+£ (ii 2 +c*-|-w 2 )J+| PaWii* 



Figube 11 .— Surfaces used in study of momentum. 

The end faces of the cylinder may be denoted, as in the 
figure, by I, II, and the curved surface by III. Then in 
supersonic flow, if a distribution of sources is restricted to a 
region downstream of I, the drag D on the body corresponding 
to the source distribution is given by the expression. 

D=jpo J u j" [(M 0 S — l)u 2 -\-o i -\-vf] dydz— poj^ j uv T rdddz 

(38) 

where v r is the radial component of the perturbation velocity. 
No loss in generality results, moreover, if the surface II is 
moved infinitely distant downstream and the radius of the 
cylinder is made arbitrarily large. The notations II and 
III will henceforth refer to this particular configuration. 
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If the drag of a lifting surface is to be calculated, the surface 
and its vortieity wake are replaced by doublet distributions 
and in that case the integral over region II in equation (38) is 
called the vortex drag of the body while region III yields the 
wave drag. (In subsonic theory, the vortex drag is custo- 
marily referred to as induced drag.) It has also been shown 
(see, for example, reference 17) that the vortex drag of a 
supersonic wing is a function only of its span load distribution 
and is equal to the induced drag at subsonic speeds for the 
same span lading. If a finite nonlifting body is considered, 
each of the velocity components in region II is attenuated 
in such a manner that its contribution to the vortex drag is 
zero. The integration over region III again gives the wave 
drag for the nonlifting body. 

The combination of the results given in this and the last 
section provides a method for finding the wave drag of an 
arbitrary body. The first step is the determination of the 
source-sink or doublet distribution corresponding to the 
body and then, by means of the principle of equivalent 
positions, the sources or doublets are moved to the x axis. 
The wave drag is then calculated from equation (38) once 
the induced velocities on the control surface are kn own. In 
the next section the wave drag will be written in a different 
form and the drag at sonic speeds will be investigated. This 
analysis also provides some insight into the range of validity 
of the sonic theory. 


V 



Figure 12. — System of axes in transformation equation (39). 


Evaluation, of wave drag as M 0 approaches one. — In order 
to study the drag of a symmetrical body at zero angle of 
attack, it is convenient to consider the general expression 
for the velocity potential given in equation (29) . Introducing 
first the transformation (fig. 12) 


where 


f=— Vi tan p+xA 
n=yi sec p ) 

tan p=f3 cos 8 


equation (29) becomes 


(39) 


4>(x,y,z)=-±j 'j 


COS p AlOp (l;,rj)d%dri 

Vte— *i) ! — P 2 (y— yiY— 


Since, however, it has been shown that <p evaluated infinitely 
far away from the wing does not change if a source is moved 
along the line |=constant, it follows that the source strengths 
can be integrated along these lines. The second integration 
is then along 17=0 where, from equation (39), f=ri and the 
value of the potential at an infinite distance is 


v, *)— ’>*> 


Setting 


f(xi, p) = — COS pf&WQ&i, T])d-q 

it follows that 


r. 


1 rx-fir 


f(x u p)dx l 


0 VCs-Zi)*-^ 


(40) 


and this is the same as the potential for a body of revolution 
with source strength per unit length given by f(xi) . The 
induced velocities corresponding to the potential in equation 
(40) are found to be, after first integrating by parts and 
using the notation p)—f'(x u p) together with the 

relation /(0, p)=0, 


b0 = J L f(x 1 , p)dxi 

dx 2irJo V(z— Zi)*— /3*r* 

. 1 C x -* T (a— zQ/'fa, p)dx l 

r= dr 2 irrJo V (x — ar r ) 1 — / 3 ®r* 


(41) 

(42) 


Asymptotic values of the velocity components for large 
values of r are readily seen to be, after first setting r=x 0 +)3r, 


u 1 C a f(x 1 , tidxi 

2rr^2f3rJo tJx 0 —Xi 


(43) 


and 


J_ IJ_ C a f{x 1 , p)dx y 
2ir V 2 r Jo ^x 0 —Xi 


(44) 


Equations (43) and (44) may be used together with equa- 
tion (38) to give for the value of drag the expression 


23 - a fis r dx, fgbiMi r 

87rJo Jo Jo \X Q —Xi Jo iJxq — Xi 


(45) 


Assuming that the body is of finite length so that _f'(x) — 0 
for reversal of the order of integration yields the relation 

D= J 0 f (x u p)f (a- 2 , p)ln\xi— x^dxrdxi (46) 


If equation (46) had been derived for a body of revolution, 
then f(x) would have been independent of the angle p and 
in that case the expression for drag would reduce to the form 

D= — f (x{)¥ (rj) ln\xi—xj[dxidxi (46a) 


This expression was given by von Karmin in reference 15. 


206 


EE POET 956 — NATIONAL ADVISORY COMMITTEE EOE AERONAUTICS 


For the study of the drag of a lifting surface, consider now 
the general expression for the velocity potential given by 
equation (30). The doublet distribution occupies in this 
case both the wing plan form and the wake since the jump 
in <j) exists also in the vortex wake. By use of the transfor- 
mations in equation (39), equation (30) becomes 


reversal of integration in equation (51) yields the simpler 
expression 

D= &0*Jo sin20 *£ X g "^ Xi> filnlxi—vildxidxt 

(52) 


6(x, y, z)=4A \ f f cos ft AjKg, y)dt-dy 

2ir IJJ [(x-xy-p '-(y-y^-pzY* 

and, exactly as in the case of the source distribution, this 
can be reduced to 


V, z) — 2 X J i, y) dy 


Setting 


it follows that 


ff(x u p) =cos p f A<j>(xi, y) dy 


r,ff) 


/3 2 r sin 0 

g(x u p) dxi 

2k 

Jo [(a— Zi) 2 — /3 a r 2 ] 3/2 


Integrating by parts and using the fact that gr(0,ft)=0, 


and 


b<f> sin.6 fx-fr ^ 2 r i g'(x 1 ,p) dx L 
U ~$x~2t r Jo [(x- Xl )*-/3 2 r 2 ] 3 ^ {47) 


t _ sin 0 5 _1 p-» (x—Xi)g'(x h p) dxi 
r dr 2irr \ r Jo /(z-a ;0 , -j8 , r* 


C «-» (x—x l )p 2 rg'(x 1 ,p)dx 1 
[Jo [(*— Xi) 2 — jS 2 r 2 ] 3/2 


(48) 


where g' (x u p) indicates ^g(x u p). Setting x =z 0 + /Sr and 

letting r approach infinity, the asymptotic expressions for 
equations (47) and (48) become, if gr'(0, ft) =0, 


u- 


-sin 6 /£ p° g"(x u p)dx l 

2ir V 2r Jo 


and 


»r = - 


sin 0 


l po g"(x i, ft) dxL 

2k ^ y 2r Jo ^J Xo — Xi 


(50) 


The relations just derived may be used in conjunction 
with equation (3'8) to give the wave drag of a lifting surface. 
This result takes the form 

»-£? r sw e r ■ ix, r c&m* p c<*. *>** (5I) 

8ir Jo Jo Jo V*o— *1 Jo ^Xq—X 2 

In the wake of aliftingwing the function p" (x u p)= 0 and if, 
moreover, 

pen - .. 

. J 0 9" fa, v)dxi=g' (c Q , p)—g'(Q, #t) =0 


It is possible to draw some general conclusions from equa- 
tions (46) and (52) regarding the wave drag of wings and 
bodies of revolution without the necessity of detailed applica- 
tions to particular configurations. It is apparent imme- 
mediately from equation (46a) that the wave drag of a body 
of revolution at zero angle of attack is independent of Mach 
number. This conclusion does not apply, however, to the 
nonlifting wing since the distribution function / (a*, p) in 
equation (46) contains the variable p which, in turn, is a 
function of both 0 and )S. As M 0 approaches one, the study 
of the nonlifting wing is divided most conveniently into two 
parts, depending on the behavior of / ( x , p). 

Consider first the more general situation in which/ (x, p) 
is not zero; that is, the case in which the number of sources 
does not-equal the number of sinks along the line £= constant. 
This means, when M 0 is 1, that an unequal number of 
sourceaand sinks appear in the transverse or yz plane and, 
if equation (46) is applied, either a finite or an infinite value 
of drag can result. The limiting value of drag at sonic 
speed, obtained from integrations of surface pressures, was 
given by Stewart and Puckett in reference (16) for several 
wing plan forms, all of which had nonvanishing values of 
/ (x, p). If the pressure distribution is calculated, however, 
the local pressure coefficients are seen to become infinitely 
large as sonic speed is reached, even for the body of revolu- 
tion, so that the assumptions of the linear theory are violated 
and the reliability of the drag predicted by equation (46) 
can in no case be assessed even though the predicted values 
remain finite. Equation (43) shows also that when control- 
surface, methods are used to compute drag at il/ 0 =l, the x 
component of induced velocity increases indefinitely when 
/ (x, ft) is not zero and that the theory is, therefore, no longer 
consistent. 

In the very special second case, that is, when / (x, p) 
vanishes for all values of 0, the analysis just presented breaks 
down at equation (40). It is clear, however, that in this ' 
case there are equal numbers of sources and sinks in the £ ±= 
constant plane and the behavior of the flow field at infinity is, 
therefore exactly the same as that which would havo been 
produced by a distribution of doublets. Equations (49) and 
(50) give the velocities induced at infinity by an arbitrary 
doublet-distribution. These induced velocity components 
are, in terms of 8, one degree higher than the similar com- 
ponents for the nonlifting case. The values of both u and v r 
can thus be expected to approach zero for all values of il/ 0 
as r approaches infinity for any flow field generated entirely 
by doublets or by an equal number of sources and sinks. It 
follows, then that the linearized theory for lifting surfaces 
(generated entirely by doublets) and for bodies with thick- 
ness distributions such that/ (x,p) vanishes (generated by an 
equal number of sources and sinks in all £= constant planes) 
is entirely consistent as M 0 approaches one and, in particular, 
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for Af* equal to one. This being true, it follows immediately 
from equation (52) that the wave drag of a lifting system is 
zero at sonic speed . 

Thickness solutions at M 0 = 1. — A swept-back wing of con- 
stant chord and infinite aspect ratio is an example of a prac- 
tical aerodynamic shape for which an equal number of sources 
and sinks occur in every yz plane. (See fig. 13.) Consider 
the case in which the wing cross section is diamond shaped 
with a slope equal to X in a plane normal to the leading edge. 
Then, in a transverse plane (section BB of fig. 13), w D equals 
± T(,X cos xp, the minus and plus signs applying, respectively, 
to the left and right of the ridge line. Accordingly, the solu- 
tion of the problem can be written in terms of a distribution 
of sources, thus 

1 f'xcoXii 

4>=—2^ f. c . T 7 oX cos P InKy-yiY+z 2 ] 



l f( I+ 2 eos i) eot * 

h+7T YoX cos $ InKy—yiy+z^dy! 153) 

■sir Jx cot <!> 


The value of bpfbx can immediately be found to be 


btj) / l r 0 Xcos^ \ 
\ 2tr tan xp ) 


Xln 


fl>H 

( x] C 0 

\ 1 1 


y-\ 

( y Co X 

1 1 1 

’+*} 

v ^2 cos xp j 

/ t a n xp_ 

\ 2 cos xp) 

tan 

j 

( 

y-{ 

X \ 
.tan xpjj 

M i! 

y - 1 

(ab). 

' 3 ) 
+Z 1 } 


(54) 


from which it is apparent that as r= z s becomes infinitely large, bcpfbx approaches zero, 
that is, for z=0, bp/bx becomes 


— I'qX cos p , 
K bxJ~ v tan p ln 


( y tan xp— x~ c 0 sec (y tan p—x+^ 

(yt&mp—xy 


Cn sec 


In the plane of the airfoil 

01 


and, using the definition for pressure coefficient, C p =—^j-, 

V o 

this becomes ' 


2X cos xp 
it tan xp 


In 



1 

r) C 0 


y sin xp— x cos 



(55) 


Equation (55) can be derived by entirely different methods. 
Perhaps the most direct of these alternative derivations is 
the one introduced by R. T. Jones in reference 18. The gen- 
eral statement used in that report is that the component of 
translational velocity of a cylindrical body in the direction 
of its long axis has no effect on the motion of a frictionless 
fluid. Hence, the pressures over the wing shown in figure 
13 are the same as those over a wing moving normal to a 
free stream with a velocity T 0 cos xp. Using the Prandtl- 
Glauert collection to the thin airfoil solution of a two- 
dimensional, diamond-shaped, nonlifting section exposed to 
a free stream with velocity T T 0 cos xp, one obtains, for 
Mo cos ^ <1, 


In l(x'-x 1 ') 2 + OS'*) 2 ] dr/ (56) 

where w 0 is the vertical induced velocity on the upper side 
of the 2=0 plane, x' is measured normal to the leading 
edge, and jS'= /l— J/ q 2 cos- xp. If this solution is referred 
to the axial system of figure 13 by the transformation 

x' =x cos xp—y sin xp 


and the integration is performed after taking the partial 
derivative with respect to x, the resultant expression for 
pressure coefficient is 


n — cos 

7T /l — Mo' 


>S 2 xp ( 2 e ' n ^ 

0 3 cos 2 xp \_\xcosxp-ysmxpJ 


(57) 


At sonic speed this equation reduces immediately to 

n 2X cos x p ( 2 Ca ) 1 

v xtan^ [_ \a; cos — t/ sin J 


(58) 


which is identical to equation (55). The result expressed by 
the two equations is, of course, not new. The significant 
point is that the same variation in pressure coefficient was 
obtained by two widely different avenues of approach and 
that the result obtained from the particular methods appli- 
cable to sonic speed theory is in agreement with that derived 
from more conventional analysis. 

lifting-surface solutions at M 0 = 1. — It should be men- 
tioned at this point that Robinson and Young (reference 19) 
have shown by means of linearized theory that supersonic 
triangular wings and subsonic elliptical wings of the same 
aspect ratio have values of lift-curve slope which approach 
a common and finite limit as M=l. The present section 
of this report is concerned only with the study of lifting 
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surfaces at a fixed sonic velocity but the results to be obtained 
are in agreement with the limiting values of reference 19. 

A further application of the results in this section can be 
made to the case of very low aspect ratio wings' at arbitrary 
Mach numbers. This viewpoint of the theory was first 
presented by R. T. Jones in reference 20 ’ and applied to 
triangular . wings, while ip reference 21^ extension was made 
to include pointed wings bn slender__bodies of revolution. 
This duality of interpretation, that is, to all aspect ratios 
at sonic speed or low aspect ratios at all Mach numbers, 
applies to all solutions of three-dimensional problems ob- 
tained from equation (13). In the subsequent analysis, 
attention will be confined to swept-back plan forms of lifting 
surfaces with pointed vertices and thus doublets will be used 
exclusively. 

In application, the two types of boundary conditions to 
be considered are as follows: 

1. Boundary-value problem, of the first kind , loading speci- 
fied . — It is given that Auo=u u —Ui = 0 over the xy plane 
except for the region occupied by the wing where 2u u = — 
2 Ui=Au 0 =f(x,y), the function being determined by the 
specified loading. Over all of the xy plane, the imposed 
conditions are Air 0 — 0. 

2. Boundary-value problem of the second kind, surface speci- 
fied . — Over the xy plane, the imposed conditions are Aw a =0 
everywhere and, except for the region occupied by the wing, 
Au 0 = 0. Over the region occupied by the wing w Q — u> u ~ 
u>i—f(x,y) where f(x,y) is determined by known camber, 
twist, and angle of incidence. (The delta notation again 
indicates the jump in the value of the variable at the 2=0 
plane. Subscripts u and l indicate conditions on the upper 
and lower surface, respectively, of this plane.) 


A4>fx, y) = ( Au 0 (x u y)dx l (60) 

J L.E. 

where the integration extends from the leading edge to an 
_ arbitrary point x. Since, moreover, load coefficient A pjq is 
related_to A u a by means of the equation 

"A p 2Au 0 

~ t ~ f 7 o 

it follows that the velocity potential <j> can bo found for any 
prescribed load distribution of a given plan form. The 
value of vertical induced velocity, evaluated at 2 = 0 , then 
suffices to calculate the twist and angle of attack of the wing. 

If the boundary-value problem is one of the second kind, 
the vertical induced velocity is given on the wing and the load 
distribution is to be found. In this case the use of equation 
(59) leads to the consideration of an integral equation. Since 
however, this integral equation is a common one in aei-o- 
dynamic theory, certain established methods may be applied 
to it. 

After noting that A<j> a (x, y)—0 at the leading edge, integra- 
tion by parts and introduction of the relation ... 


yields for perturbation potential the expression 

Ai> 0 (x, yi) arc tan V ~J y dy x (61) 

In the limit as z approaches zero the derivative of 4> with 
respect to z reduces to the form 


y 



The nature of the differential equation shows that the 
value of <j> is a consequence of boundary conditions along 
lateral strips. If, as in figure 14, the two leading edges are 
given by the expressions 1 y — bi {x) and y = b 2 (x), the velocity 
potential is expressible in the form .... 


4>(z, y, 


, _z_ pa A<j> 0 (x, y^dyi 

' 2ir Jb, {y—yiY+z* 


(59) 


If the boundary-value problem is one of the first kind, the 
general expression for ® follows from a direct integration 
after noting that : - ~ - - 


1 f». ' (62 ) 

2 it J ft, y~y x _ 

For a given distribution of w 0 over the plan form of the 
wing, equation (62) represents an integral equation to be 
solved for Av 0 (x,y) subject only to the condition that the 
Kutta-Joukowski condition is satisfied at all subsonic trailing 
edges. Once Avfx,y) is determined it follows that 


A0o=J^ Av^x^dyx 

(63) 

A p 2 5A<ft) 

q Vo 

(64) 


In the present report the solution to the wing plan form 
shown in figure 15 will be presented^ The value of A<j> a 
which satisfies equation (62) is, in region 1, 

A0o=— 2w 0 fx s tan* B—y 2 (65) 

and, in region 2, __ 

A<j> 0 =—2wyx tan 8 [E(\p 0 ,k 0 ) —k 0 ' s F (tafia)] (66) 
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where E and F are defined in the appendix and where 

, . j x 2 tan 2 6—y 2 

f,=.ro ^ 

dl 


k 0 r =‘-^— = ^IT^h i 

x tan 6 


C67) 


( 68 ) 


The equation y—ai{x) of the trailing edge for which equation 
(66) is valid is given by the formula 


at— 


k ' 

E 0 -h n K a 


(69) 


which expresses a L explicitly as a function of ( - - - Sl ~ - )- 

\2? t- 8iB u J 


Regions 



This particular choice of trailing-edge shape was used to 
simplify the analysis. The resulting plan form approaches 
a constant-chord wing as the span increases. The variation 
of di with x is given in figure 16; and figure 17 shows the 
relation between aspect ratio and span. 

The loading coefficient is given in the two regions (defined 
in fig. 15) as follows : 




Region 1. 


Region 2. 


A p 4z tan 2 8 

qa Vz* tan 2 6—y 2 


(70) 


(71) 


This load distribution is shown in figure 18 for a triangular 
and a swppt-back wing. It is seen that the loading at sonic 
speed bears a close resemblance to those found at higher 
Mach numbers. Two similarities of note are the discon- 
tinuity in the pressure gradient at the Mach wave originating 
from the trailing edge of the root chord and the satisfying 
of the Kutta condition only where the trailing edge is sub- 
sonic. The lift and vortex drag coefficients are given, 
respectively, by 



Figure IS.— Pressure distributions for triangular and swept-baek wings at Me=l. 


and 


Op, _ A py- EZ-k^Ks' l 

a? tan 6 tan 0 [_ 4 (s 0 / c o tan 0) J 


( 73 ) 




REPORT _956 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


210 

where 

V=J=VI=£7 

These coefficients are plotted as a function of Ajt&n 6 in 
figure 19. It is shown that the values of C D Jo? tan 8 and 
Cz.la tan 6 for finite aspect ratio swept-baclc wings are always 
less than the corresponding values for the triangular wing 
(A/tan 0=4). When the span of the swept-baclc wing 
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Figure 16.— Variation of lift and drag with aspect ratio for a swept-back wing at Mq=1. 

i 

becomes very large, the slope of the trailing edge approaches 
asymptotically the slope of the leading edge. It follows 
that for infinitely large aspect ratio the limiting value of the ' 
load distribution on the outboard sections should approach 
the value given by simple sweep theory for an infinitely long 
swept-back lifting surface with constant chord. .This result 
is, in fact, a consequence of equation (71). 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffet Field, Calif., Dec. 12, 194-8. 
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